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Abstract 

We deal with the scattering of an acoustic medium modeled by an index of refraction 
n varying in a bounded region Q of and equal to unity outside 0. This region is 
perforated with an extremely large number of small holes Dm’s of maximum radius a, 
a « 1, modeled by surface impedance functions. Precisely, we are in the regime described 
by the number of holes of the order M := the minimum distance between 

the holes is d ~ a* and the surface impedance functions of the form Am ~ Am,o«~^ 
with /3 > 0 and Am,o being constants and eventually complex numbers. Under some 
natural conditions on the parameters /3, t and Am, 0 ; we characterize the equivalent medium 
generating, approximately, the same scattered waves as the original perforated acoustic 
medium. We give an explicit error estimate between the scattered waves generated by 
the perforated medium and the equivalent one respectively, as a —)■ 0. As applications of 
these results, we discuss the following hndings: 

1. If we choose negative valued imaginary surface impedance functions, attached 
to each surface of the holes, then the equivalent medium behaves as a passive 
acoustic medium only if it is an acoustic metamaterial with index of refraction 
h{x) = —n(x), X € Q and n(x) = 1, x G \ fl. This means that, with this 
process, we can switch the sign of the index of the refraction from positive to nega¬ 
tive values. 

2. We can choose the surface impedance functions attached to each surface of the holes 
so that the equivalent index of refraction fi is n{x) = 1, x G M^. This means that 
the region modeled by the original index of refraction n is approximately cloaked. 
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1 Introduction 

The derivation of the macroscopic behavior of a given physical system (or other biological, 
chemical systems, etc.) as an interaction of a ’dense’ microscopic particles is a well known 
procedure for a long time, see mu- However, the mathematical modeling as well as the 
justihcation of this procedure was understood only in the middle of the last century, see [TJ 
ESI ED]. One of the most known ideas to describe the passage from the microscopic states 
to the equivalent macroscopic state is the homogenization theory, see [H ESI ED]. There are 
two approaches: the deterministic one and the probabilistic one. If one wants to estimate the 
macroscopic state deterministically, then one needs to assume the periodicity in distributing 
the small particles. To avoid the periodicity, one can assume the small particles to be randomly 
distributed and in this case we estimate the equivalent medium in the probabilistic sense. Let 
us emphasize here that in both approaches, we estimate the limit of the helds created by the 
microscopic structure to the helds created by the macroscopic one with energy norms stated in 
the whole domain where the small particles are distributed. 

If we are interested only in estimating the limit of the helds away from these small particles 
(as in the inverse problems and the design theories), then another alternative, avoiding both 
the periodicity and the randomness, is possible. 

A root of this alternative goes back to the seminal work by L. Foldy, see 021 , where he gave 
a close form of the held scattered by M isotropic point-like scatterers, see [HIEII for details on 
this model. The justihcation, or the mathematical foundation, of Foldy’s representation was 
proposed by Berezin and Faddeev in another seminal work, see |D|. The idea is that based on 
the Krein extension theory of self-adjoint operators, one can model the dihusion by point-like 
particles by the Schroedinger model with singular potentials of Dirac type supported on those 
point-like scatterers. This has opened a very fruitful direction of research, see the book |2|, on 
exact models. Following Fadeev’s approach, using Krein’s extension theory, the dihusion by 
small particles was studied in several works, see for instance [IHl ESj. 

A diherent, but still related, approach to describe the dihusion by small particles is based on 
the integral equations. This is suggested by diherent authors, as A. Ramm 1271 and H. Ammari 
and H. Kang [I], for instance. In particular, A. Ramm [271 EH] shows that the dominant term 
in the expansion of the scattered held has the same form as the Foldy’s close form where the 
centers of the small scatterers play the role of the point-like particles. He used the (rough) 
condition | ^ 1, where a is the maximum radius of the small particles and d the minimum 
distance between them, and no error estimate is derived, but he formally characterized the 
equivalent medium. We also cite the recent works by V. Maz’ya and A. Movchan [22l ES] where 
they study the Poisson problem and obtain error estimates. In their analysis, they rely on the 
maximum principle to extend the boundary estimates, which argument does not go smoothly 
for stationary models as the Helmholtz one. 

A rigorous approximation of the scattered fields, with error estimates in terms of the number 
M of the small scatterers, their minimum distance d and the maximum radius a, was derived in 
[DiEni using the integral equations approach. Based on these last estimates, we can characterize 
the equivalent medium with explicit error estimates in terms of the parameter a, a « 1, 
in appropriate regimes described by the other parameters M and d in terms of a, namely 
M := M(a) := 0{a~^) and d := d{a) ~ a*, as a « 1, with non negative parameters s and t. 
This was done in where the small particles are taken to be soft acoustic or rigid elastic 
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particles. The objective of the present work is to extend this stndy to the case where the small 
particles have impedance type snrfaces with the scaled snrface impedences of the form ~ 
with non negative /?, in addition to the scaled coefficients M and d described above. Compared 
to the works in [1], we derive here better error estimates. Precisely, fixing s = 1, /d = 0 and 
f = I for simplicity and as an example, we derive here an error of the form 0 ( 0 ™“'^'’'’ sf) while 
in [1] it is of the form Here 7 G (0,1] is the Holder regnlarity exponent of the 

coefficients Aq and K appearing in the equivalent medium, as discussed below. 

The design of materials with desired, and in particular negative, index of refraction is a 
hot topic in the last years, see [H] for instance. Concerning this topic, our contribution is to 
have shown mathematically that this is possible with high generality. Indeed, we show that 
the equivalent medium is modeled by three coefficients K, Pq and Aq modeling respectively, the 
local distribution of the particles, their geometry and the impedance coefficient attached to 
each particle. Since we have the freedom in choosing the three functions K, Pq and Aq, then we 
can generate a large family of indices of refraction. In particular 

1. if we choose the surface impedance to have negative imaginary parts, which is mathemat¬ 
ically possible, see [lO], then we show that the equivalent medium will be passive only if 
the index of refraction is negative, i.e. it behaves as a metamaterial. 

2. we can choose the coefficients K, Pq and Aq so that the equivalent index of refraction is 
educed to the unity. This means that the domain H modeled by the original index of 
refraction n is cloaked. 


Let us emphasize that the derived explicit error estimates between the fields generated by 
the microscopic structure and the one generated by the equivalent macroscopic structure might 
be useful to quantify the accuracy of the design. 


Let us also cite some related works on the derivation of the equivalent media. The first works 
go at least to Rauch-Taylor, see |29j, see also the works by Cioranescu and Murat [Ellin], who 
characterized the limiting problems for some Poisson type problems and provided convergence 
results (but with no rates of convergence) of the corresponding resolvent operators. Later, 
these results were extended and refined in the works of Ozawa and the ones of Figari et ah, see 
respectively |2B] and [TH] for instance, using point interaction approximations of the Green’s 
kernels. Compared to these results, we do not need the periodicity nor the randomness in 
distributing the small particles in addition we model them via the scaled parameters M, d and 
the Am’s with a high generality as described above. 


The rest of the paper is organized as follows. In section we state the main results. 
Precisely, we describe the mathematical model of the stationary scattering by many small 
bodies of impedance type in section 2.1, then we state the main mathematical results in section 
2.2 We end this section with a discussion on the possible applications of these results in the 


acoustic metamaterials and the acoustic cloaking in section 2^ Section is devoted to the 
proof of the main theorem of the paper. 
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2 Statement of the results 


2.1 The acoustic scattering by many impedance type holes 

Let i?i, i? 2 , • • •, Bm be M open, bounded and simply connected sets in with Lipschitz bound¬ 
aries containing the origin. We assume that the Lipschitz constants of Bj, j = are 

uniformly bounded. We set Dm ■= eBm + Zm to be the small bodies characterized by the pa¬ 
rameter e > 0 and the locations Zm fn = 1,, M. Let 17* be a solution of the Helmholtz 

equation (A -|- k‘^)W = 0 in We denote by the acoustic held scattered by the M small 
bodies Dm C due to the incident held f/* (mainly the plane incident waves W{x, 6) := 
with the incident direction 6^ G where being the unit sphere), with impedance bound¬ 
ary conditions. Hence the total held f/* := 17* -|- 17® satishes the following exterior impedance 
problem of the acoustic waves 


M - 

(A + KV(a:))17* = 0 inM"\( UD, 

V m=l 


dU^ 

dl^m 

dU^ 

d\x\ 


+ 


= 0, 1 < m < M, 


dD„ 


inll^ ~ ' I—i ) ’ 1^1 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


Again, the scattering problem (2.1 2.3) is well posed in the Holder or Sobolev spaces, see 


[II1DS1I21] in the case ‘^Xm > 0. As we said for ( |2.1||2.3[ ), this last condition can be relaxed 
to allow '^Xm to be negative, see [10]. Applying Green’s formula to 17®, we can show that the 
scattered held 17® (x, 6*) has the following asymptotic expansion: 


U^{x,9) = 


giK|x| 

47r|a; 


■U°°{x,6) + 0{\x\ ), |a;| —)■ oo. 


(2.4) 


where the function 17°° (7, 6) for (7, 0) G x is the corresponding far-held pattern. 
Definition 2.1. We define a := max diam{Dm), d := min dmj where dmj '■= 

l<m<M 

dist{Dm, Dj) and set Kmax as the upper bound of the used wave numbers, i.e. k G [0, Kmax]- 
The distribution of the scatterers is modeled as follows. 

1. The number M := M{a) := 0{a~^) < Mmaxa~^ with a given positive constant Mmax- 

2. The minimum distanced := d{a) ~ of, i.e. dminoX < d{a) < dmaxoX, with given positive 
constants dmin and dmax- 

3. The surface impedance Xm '■= Xmpa~^, where Xm,o 7^ 0 and might be a complex number. 
Here the real numbers s, t and jd are assumed to be non negative. 

We call the upper bounds of the Lipschitz character of Sm’s, Mmax, dmin, dmax and Umax fhe 
set of the apriori bounds. In [10], we have shown that there exist a positive constant oq, A_ 
and A+ depending only on the set of the apriori bounds and on Umax such that if 


a < ao, |Am,o| < A+, |3fJ(Am,o)| > A_, /3 < 1, s<2- jd, -<t 


(2.5) 
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then the far-held pattern U°°{x,6) has the following asymptotic expansion 


M 


= KTli.s) + Y1 + O (o=>—"'’), 


( 2 , 6 ) 


m=l 


uniformly in x and 0 in The constant appearing in the estimate 0{.) depends only on the 
set of the apriori bounds, A_, A+ and on Umax- The quantity —x) is the total held 

evaluated at the point Zm in the direction —x, corresponding to the scattering problem 


(A -I- K^n^{x))Vjl = 0 in 


dV^ 

d\x\ 


inV^ = o ( , |a;| —)• oo, 


\x\ 


(2.7) 

( 2 . 8 ) 


i.e. V^{zm, —x) := e -|- V^{zm, —x), where is the scattered held. The coefficients 
m = 1, are the solutions of the following linear algebraic system 


M 


CmGK,{Zm, Zj)Qj — —CmVj^{Zm,0), 


(2.9) 


i=i 


for m = where Cm ■= —Am \dDm\- Here Gf^{x,z) is the outgoing Green’s function 

corresponding to the scattering problem ( |2.7 2.8). 

The algebraic system (2.9) is invertible under the condition: 


s < 2 - /3. 


( 2 . 10 ) 


2.2 The equivalent model 


As the diameter a tends to zero the error term in (2.6) tends to zero for t and s such that 

/3<1, s<2-A (2.11) 

and it is at least of the order Observe that we have the upper bound 


M 


Ve-“Qm|<M sup \Qm\ = 0{a^-f^-n 

m=l ’ ’ 


( 2 . 12 ) 


since Qm ~ |Am||T)m| ~ see [lU]. Hence if the number of holes is M := M(a) := 

(2.13) 


0{a ^), s < 2 — /3 and t satishes (2.11), a —)■ 0, then from (2.6), we deduce that 
U°^{x, 6) —)■ V^{x, 9), as a —)■ 0, uniformly in terms of 9 and h in 


This means that this collection of holes has no effect on the homogeneous medium as a —)■ 0. The 
main concern of this work is to consider the case when s = 2 — /3. Let H be a bounded domain, 
say of unit volume, containing the holes Dm, m = 1,..., M. We divide H into [a^~^] sub-domains 
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Qm, ITT' = I,---, such that each Q^n contains D^, with Zm € flm as its center, and some 

of the other -Dj’s. We assume the number of holes in flrn, for m = 1 ,to be uniformly 
bounded in terms of m. To be precise, we introduce K : —>■ M as a positive continuous 

and bounded function. Let each him, m G N, be a cube such that Qm H hi (which we denote 
also by Qm) is of volume contains [K{zm) + 1 ] holes (where [a] stands for the 

integral part of a G M). We set Kmax ■= sup 2 ^(iL(zm) + 1), hence M = Y^j=i + 1] < 

Kmax[a^~‘^] = 0(a^“2). 



Figure 1: An example on how the holes are distributed in hi. 

We prove the following result. 

Theorem 2.2. Let Aq : hi '—> C be a continuous function and take Xm ■= Xo{zm)(i~^- Consider 
the small holes to be distributed, as described above, in a bounded domain hi, say of unit volume, 
with their number M := M(a) := and their minimum distance d := d{a) := of, 

< t < 2 — /3 with (3 < 1, as a ^ 0. 

1. If the shapes of the holes are different, under the condition that the reference bodies Bm’s 
have uniformly upper bounded perimeters with uniformly lower bounded radiis, then there 
exist a function Pq in r]p>iLP(Mf) with support in hi such that 

limt/“(x, 6 *) = Uf°{x,9) uniformly in terms of 6 and x in (2-14) 

a —)-0 

where U^{x,6) is the far-field corresponding to the scattering problem 


(A + + {K + l)PoAo)t/o = 0 in M'^, 


U^o = U^ + e 

■ rrs ( ^ 

— -%kUq = o — 
o\x\ \|a:| 


IKX-d 


\x\ —)■ oo. 


(2.15) 

(2.16) 

(2.17) 


2. Assume in addition that Aq and K |q are in C°’'^(hl), 7 G (0,1] and the reference bodies 
Bra’s have the same perimeter and diameter, and denote by P := Then 

U°°{x,e) = 17o“(L,0) + 

uniformly in terms of 0 and x inS^, where Pq — P in Q and Pq — 0 m \ fl. 


(2.18) 
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We see from case 2 of Theorem 12.21 that 


u^{x,e) = u^{x,e) + 


V- 2-/3-t}) if 0 < /3 < I 
1-3/3, 2-^-t}) if | < /? < f 


(2.19) 


and for (3 > the remainder is no longer tending to zero as a tends to zero. 

We also see that the best error estimate is attained for (3 = 0 and it is 5’ In 

particnlar if we reasonably take t < 1, which means that the minimnm distance is of the order 
of the diameters, i.e. d ^ a, and 7 > then 

U°°{x, 6) = f/o“(£, 6) + 0(ai), a ^ 0. (2.20) 


However in this case, i.e. (3 = 0, the nnmber of small holes attains its maximnm, M = 
= 0(a~‘^). Actnally, there is a compromise between the nnmber of the small holes and 
the order of the approximation. In short, the larger is the nnmber of the small holes (or the 
smaller is (3) the better is the approximation. 


2.3 Some applications 

As a corollary of Theorem m we dednce the following resnlts. 

1. We write + {K + l)Po-^o = he. This way of representing 

the eqnivalent coefficient + [K + l)Po-^o means that the eqnivalent material behaves 
as an acoustic material whose index of refraction is h satisfying 


particular, we can choose Aq, see remark 2.3 


as 


•^0 •— -^0 


( 2 . 21 ) 


then 


= n‘^ + {K + l)PoXo. (2.22) 

We set h = hi + m2. This new acoustic material will be passive only if A h = h2 > 0. 
From (2.22), we deduce that 


n. 


hi = + 2tt{K + l)Po3^ Aq and hi 77,2 = vr(iF + l)PoA Aq. 


(2.23) 


Recall that the coefficient Aq comes from the surface impedance functions Am attached 
to every small body of the collection of the small bodies generating the coefficient + 
{K + l)PoAo. Hence, if we choose Am’s so that AAq > 0, then necessarily hi > 0 and then 
we can generate acoustic materials having index of refraction as 


> 0 and An > 0. 


(2.24) 


Now, if we choose Am’s so that A(Ao) < 0, then we deduce from ( 2.23[ ) and the fact that 
h 2 > 0 that necessarily hi < 0. With this way, we can generate acoustic materials of the 
form: 


3? n < 0 and A n > 0. 


(2.25) 












In addition, if we choose the surface impedance function so that Aq := Ao(e) with 

(e)( < 0 ) —)■ 0, as e —)■ 0, so th at the cond 
, then the two equations in (2.23) imply that 


is satished i.e. f! > 


3? Ao(e) ( > 0 ) —)■ 0 and 9 Ao(e)( < 0 ) —)■ 0, as e —)■ 0, so that the condition (3.9) 

_ VMM, 

|Ao(e 


3? n{e) —>■ —n and O' n(e) —>-0 as e —)■ 0. 

As a conclusion, if we perforate a given acoustic material, modeled by the index of re¬ 
fraction n{x), with appropriately distributed small holes with negative imaginary part 
surface impedance functions, attached to each surface, the equivalent medium behaves 
as a passive acoustic medium only if it is an acoustic metamaterial, i.e. with index of 
refraction 

h{x) = —n{x). (2.26) 

2. We can choose the surface impedance functions so that =1, x G 12. For 

instance we take Aq := Aq and Aq := (x+i)Pp • Hence the equivalent index of refraction h 
is n{x) = 1, X G This means that the region modeled by the index of refraction n is 
approximately cloaked. Observe that since we can take the surface impedances complexe 
valued with any sign then we can cloak the region O dehned by complex valued index of 
refraction with any sign of the real and imaginary parts. 

Remark 2.3. The surface impedance in (2.21), i.e. Aq := Aq will be achieved if we choose 


the surface impedances of the small holes as 

Am := Am(K) := Am,o (2.27) 

with a « 1, for instance. Since these surface impedance functions appear in the boundary 


conditions (2.2), i.e. 


dUr, 

dOrr. 


+ Xm{f^)Um = 0, 1 < m < M, on dD^, 


which we can rewrite for convenience, to link the acoustic pressure Um to the velocity on the 


boundary as 


dVrr 


Um (TmiU) 


dU^ 

dVrr. 


= 0, 1 < m < M, on dDr 


(2.28) 


where cTm(«^) := A,^^(k). In the time domain these impedance boundary conditions are translated 
as 


r°° dU 

Um{t,x)+ / (Jm{t-t')-^UX!L[tfx)dt' = Q,xondDm- 

J —oo 

where d^{t) := A Jf°^a{K)e-^^*dK and Um(t,x) := A Jf°^U^{K,x)e-^'^^dK. 


(2.29) 


The wave propagation with the type of time domain impedance boundary conditions in (2.29) 


were recently object of studies see for instance m- Sufficient conditions on the admissibility of 


the impedance boundary conditions (2.29), as the reality, passivity and causality conditions, are 


given in the book fSUf and also discussed in fWj- These conditions are given in the freguency 
domain and, ifift Aq > 0 and due to the decay in terms of k, our surface impedances (JmiU) = A 


-1 

m 


with Am given in (2.21) satisfy those conditions. As a conseguence, we do hope that the choice 
of the surface impedances used in the applications described above might make sense in practice. 
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3 Proof of Theorem 12.2 


3.1 The relative distribution of the small bodies 

We start with the following observation from [3] on the relative distribution of the small bodies. 
For m = 1,... ,M hxed, we distinguish between the obstacles Dj^ j ^ m, by keeping them 
into different layers based on their distance from Dm- Let us hrst assume that K{zm) = 0 
for every Zm- Hence each Dm has the (same) volume and contains only one obstacle D 
We arrange these cubes in cuboids, in different layers such that the total cubes upto the n' 

layer consists of ( 2 n + 1 )^ cubes for n = 0 ,..., ~ l) ]> located at the center, 

see Fig 


m* 

th 


2.2 


Hence the number of obstacles, we denote by D'^ and located in the n*", n 7 ^ 0, 
layer will be [( 2 n + 1 )^ — ( 2 n — 1 )^] = 24n^ + 2 and their distance from Dm is greater than 


2-/3 

n I a 3 — ^ 

estimate 


Observe that, < (aO'' — ^) < 


2-g 


2-/3 

a 3 . 


d{D^,Dm) > 


2-li 

na 3 


Hence we deduce the needed 


(3.1) 


Now, we come back to the case where K{zm) 7 ^ 0. As | < < 1, then with such Om’s, 

the total cubes located in the layer consists of at most the double of [( 2 n -|- 1 )^ — ( 2 n — 1 )^], 
i.e. -|- 4 and the inequality (3.1) is also verihed. 



Figure 2: Rubik’s cube consisting of two layers 


3.2 Solvability of the linear-algebraic system (2.9) 


The algebraic system (2.9) can be written in compact form as 

BQ = 

where Q, U'^ G and B G are dehned as; 


B : = 


/ —^ —Gk{zi,Z2) —Gk.{zi,Z3) 

— Gfi{z2,Zi) —^ —Gk,{z2,Z3) 

\ -G^{zm,Zi) -G^{zm,Z2) 


-Gf,{zi,ZM) \ 

-G^{z2,Zm) 


Gi^(^zm: Z][/[—\^ Q J 


(3,2) 


(3,3) 
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Q:=(Qi Q 2 ... Qm and := ( ... V;‘(^m))^. (3.4) 


The following lemma provides us with the needed estimate on the invertibility of (3.2). 


Lemma 3.1. ITe distinguish the following two cases: 


> 


22-/9 


Let 3fJ(Amo) < 0 and assume that min 

l<m<M 

and the solution vector Q of (3.2) satisfies the estimate 


\Cm\^ 


\/2Mrr 


then the matrix B is invertible 


and then 


M 


M 


^ 




m=l 


/ min iR.Cm 

I l<m<M 

1 max ICmP Ti a?~h 

\ l<m<M 


-2 


M 


m=l 


(3.5) 


E l'3"*l S 2 

m=l \ 


/ min mCr, 


m max |Cm|\/26M„ 

l<m<M l<m<M 


M 


max I Cm I 

\ l<m<M 


71 a 


2-/3 


(3.6) 


m=l 


min 

• Let 3fJ(Am,o) > 0 and assume that 1 x 2 > 

T<m<M 


^‘2Mmax 

iTa?~P 


and the solution vector Q of (3.2) satisfies the estimate 


then the matrix B is invertible 


M 


EI<3™I"S4 

m=l \ 


/ min C^m) /o n/r 

/ l<m<M yllVlr, 


M 


max ICmP 77 a^~h 

\ l<m<M 


Eh« 

m=l 


(3.7) 


and then 


M 


<2 I 


min 3fJ—Cm max |Cm|\/26M„ 

l<m<M l<m<M 


M 


max I Cm I 

\ l<m<M 


71 a 


2-/3 


l<m<M ' ' 


m=l 


(3.8) 


The proof of this lemma is given in [10] for the case where n = 1, i.e. C^ is the fundamental 
solution ^i^{x,y) := of the free Helmholtz model. But that proof goes smoothly for 

as well. 

min |5RCm.| .—— 

Since Cm := — Aml^^ml, the condition 1 ^ -p > is satished if A_ and A+ 

l<m<M ^ 


satisfy 


A- ■\/26Mmax 

A+ ^ 71 


(3.9) 
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3.3 The limiting model 

From the function iF, we define a bounded function Ka : —)■ M as follows: 


Ka{x) := Ka{Zm) ■ = 


K(Zm) + l if 


X G flrn: 


0 


if X ^ flrn for any m = 1,..., [a^ ^]. 


(3.10) 


\dBm\ 


Hence each flm contains [Ka{zm)] obstacles and Kmax ■= Ka{zm)- 

Let Ca be the piecewise constant functions such that Ca\u^ = Cm '■= yi 

(max diam{Bm)] 

for all m = 1,..., M and vanishes outside fl. The constants Cm are independent of a. We set 


C := max \Cm\oo- 

l<m<M 


Consider the Lippmann-Schwinger equation 


Ua{z)+ / GCz,y)Ka{y)Ca{y)Ua{y)dy = -V^{z,e),z en 

Jn 


and set the Poisson potential 


V{Y){x) := / GCx,y)Ka{y)Ca{y)Y{y)dy, 
Jn 


X G 


(3.11) 


(3.12) 


(3.13) 


The coefficients Ka and are uniformly bounded. The next lemma concerns the mapping 
properties of the Poisson potential. 

Lemma 3.2. The operator V : —)■ is well defined and it is a hounded operator 

for any hounded domain Q in i.e. there exists a positive constant Cq such that 


l|r(r)l|a.,„, < C(,P1U>(n). 


(3.14) 


‘i\x — y\ 


Its proof is given in|T2], for instance, in the case when G^x^y) = ^KC,y) := 4 ^| 3 ,_y|- 
However satishes (A + — $«) = k^(1 — n‘^)^K, in By interior estimates, we 

deduce that ||Gk(-, z) — <!)«,(•, z)\\H2{n) is uniformly bounded in terms of z. By this last property, 
we show that Lemma 3.2 is also valid for n ^ 1 and n = 1 in \ f2. 


Using Lemma 3.2, the fact that the operator I + V : Lfi{VL) \—y Lfi{VL) is Fredholm with 
zero index and the uniqueness of the scattering problem corresponding to the model 


(A + kV - iFaCa)U = 0, in 


(3.15) 


(where Y := + and Y^ satishes the Sommerfeld radiation conditions and U* is an incident 

held), we have the following lemma, see [1] the details. 


Lemma 3.3. There exists one and only one solution Y of the Lippmann-Schwinger equation 

(3.16) 


(3.12) and it satisfies the estimate 

l|r|lL»(n) < C||K;||„.,n, and ||Vy|U»|„) < C'||l/'ll ir2(Q), 


where being a large hounded domain which contains fl. 
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3.3.1 Case when the obstacles are arbitrarily distributed 

From the definition of C^, we have €„ Ao(^m)• AVe have assumed the 

reference bodies -Bm’s to have different but uniformly bounded perimeters with uniformly lower 
bounded radiis. In addition, since Aq is a continuous function, the function Xa : —> C 
defined A^ converges to Aq uniformly. Then there exists a function Pq in L^(Q) 

such that Ca converges weakly to Cq := AqPo in L^(Q). Now, since K is continuous hence Ka 
converges to {K + 1) in and hence in Then we can show that KaCa converges 

to {K + l)Co in 

Since KCa is bounded in then from the invertibility of the Lippmann-Schwinger equa¬ 

tion and the mapping properties of the Poisson potential, see Lemma |3.3[ we deduce that 
||h^alli^ 2 (o) is bounded and in particular, up to a sub-sequence, t/* tends to Uq in From 

the convergence of KaCa fo {K -|- l)Co and the one of 17* to Uq and 
following equation satisfied by Ul{x) 

Ul{x)+ [ G^,{x,y){Ka){y)Co{y)UQ{y)dy = -V^{x,e) in fl. 

Jn 

This is the Lippmann-Schwinger equation corresponding to the scattering problem {A + — 

{K + 1)Co)17q = 0 in Uq = + U\ and satisfies the Sommerfeld radiation conditions. 

As the corresponding far-fields are of the form 

U^{x,e) = f G^{x,y){K + l){y)Co{y)U^,{y)dy 
Jn 

which can be written, by the mixed reciprocity relation <7“(7, y) = Vn{y, —x), as 

U^{x,e)= [ Vr.{y,-x){K + l){y)CoiyMiy)dy 

Jn 

and similarly the ones of t/* are of the form 

UT{x,e)= [ Vn{y,-x)Ka{y)Ca{y)U^Mdy 

Jn 

we deduce that 

U^{x, 9) — U^{x, 9) = o(l), a —)■ 0, uniformly in terms oi x,9 e 

3.3.2 Case when K is Holder continuous 

If we assume that K G 7 G (0,1], then we have the estimate ||(iF + 1) — Ka\\L°°{n) A 

Ga^, a « 1. Similarly, since we assume the shapes to have the same perimeter and diameter 
and Ao to be in 7 G (0,1], then we have also ||Co — Ca||Loo(Q) < Ga^ since Co = 

Ao ''^hsre B is the common reference body. Since the obstacles have the same perimeter, 

we set Co to be a constant in and Co = 0 in \ 12. Recall that Uq and Ua are solutions of 
the Lippmann-Schwinger equations 

Uo+ [ G,{x,y){K + l){y)Co{yMiy)dy = v:; 

Jn 


(3.12), we derive the 
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and 


Ua+ / G,{x,y)KM<^,{y)Ui{y)dy = Vl 
Jn 

From the estimate \\{K + 1) — Ka\\L°°{Q) < Ca^, a « 1, we derive the estimate 
U^{x, 6) — U^{x, 6) = 0{d^), a << 1, uniformly in terms of x, 6 e 


(3.17) 


3.4 The approximation by the algebraic system 


For each m = 1,..., M, we rewrite the equation (3.12) as follows 

M 

Ua{Zm) + G^{Zm, Zj)CjUgjzj)a^~^ 


i=i 

j^m 


M 


-V^iZrn, 0) + J2 ^j)CjUa{z,)a^-f^ - Y, G.{Zm, Z,)Ka{z,)C,Ua{z,)Vol{n,) 


i=i 


1=1 

j^m 


Y G^{z^,Zj)Ka{zj)GjUa{zj)Vol{Qj) - / G^{z^,y)Ka{y)Ca{y)Ua{y)dy.{3.18) 


1=1 

iT^m 


Let us estimate the following quantities: 


7./3-21 


4:= G^{zm,z,)Ka{zj)C,Ua{zj)Vol{nj) - / G^{z^,y)K^{y)Ca{y)U^{y)dy 


and 


1=1 

l^m 


M 


B ■=Y^'^i^rn,Zj)GjUa{Zj)a^ ^ - ^ G ^,{Zm, Zj)Ka{Zj)GjUa{Zj)V ol{VLj). 


1=1 

l^m 


1=1 

l^m 


3.4.1 Estimate of A 


/3-21 


By the decomposition of fl, fl := fl;, we have 


7 / 3-21 


G^{zm,y)Ka{y)Ca{y)Ua{y)dy = Y / G^{z,n,y)Ka{y)Ca{y)Ua{y)dy. (3.19) 

i=i 


Hence, ^■= G^{zm,y)Ka{y)Ca{y)Ua{y)dy 

r^^-2i 


E 

1=1 

l^m 


G^{z^, Zj)Ka{zj)GjUa{zj)Vol{Qj) - / G^{zm,y)Ka{y)Ca{y)Ua{y)dy 


(3.20) 
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For I ^ m, we have 

[ y)Ka{y)C,{y)Ua{y)dy - G,{z^, zi)Ka{zi)CiU^{zi)Vol{Qi) 


d^a(^^i)Gi I \Gi^(^ZjYi, y^Uai^y^ z ;)? 7 q ,(^ z )] 

J 0,1 

(3.21) 

We set f{zm,y) = G^{z„T,,y)Ua{y) then f{zm,y) satishes 

fiZm, y) - fiZm, Zi) = {y- Zi)R\Zm, y) 

where 

R\zm,y) = [ Vyf{zm,y - Piy - zi))dp 

Jo 

= [ '^y[GK{zm,y - P{y - zi))Ua{y - P{y - Zi))] dp 

Jo 

= [ [^yGnizn^y - P{y - zi))]Ua{y - p{y - Zi))dp 


+ / G^{Zm, y- P{y- Zi)) \SyUa{y - P{y- Zl))] dp. 


(3.22) 


We set := ^ iT,\x-y \ —• have the following lemma: 

Lemma 3.4. We have the asymptotic expansion: 

d^.G^{zj, y) = y) + 0(1), for y in n, with a = 0,1. 

iKn(zj)\x-y\ 


(3.23) 


Proof We know that (A + Hpn‘^{x))GK = —5{- — y), in and ^K,j{x, y) := iT,\x-y \ — satishes 
(A + n?n^{zj))^K = —dp — y), in where both G^ and <Fkj satisfy the Sommerfeld radiation 
condition. Then Hf^j{x,z) := (G^ — ^fij){x,z) satishes the same radiation condition and 

(A + Hpn‘^{zj))Hf^j = rppnp{zj) — n^{x))GK,, in M^. (3.24) 


Multiplying both sides of (3.24) by <Fkj and integrating over B DD Vt we obtain: 

H^j{x,y) = k‘^ {n‘^{zj) - n‘^{t))^^j{t,x)G^j{t,y)dt 

Jb 


(3.25) 


+ / H^j{t,y)dy{t)^f^j{t,x)ds{t) + / di,{t)Hnj{t,x)^,,j{t,x)ds{t), 


IdB 


'an 


and then 


V:,H^j{x,y) = {n‘^{zj)-n‘^{t))G^{t,y)Va;^^j{t,x)dt 


(3.26) 


IB 


IdB 


H^j{t,y)Vxdu{t)^K,jit,x)ds{t) + / du{t)H,,j{t, y)Vx^njit, x)ds{t). 


IdB 


Since {n‘^{zj) — n‘^(t))'Vx^K(t, Zj) = 0{\t — Zj\),t G B, as the singularity of Vx^ndityX) is of 
the order \t — x\~‘^, then both integrals appearing in (3.25) and (3.26) are of the order 0(1) for 
X = Zj and y E fl CC B. □ 
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From 

the explicit form of 

$ 

^ K^nri’) 

we 

lave 

^K,m{x,y) 


^ V- 

Now from 

Section 

3.1 


we Vy^K,ra{x,y) = 

precisely the inequality 


(3.1), we see that for / 7 ^ m 


c c 

\^K,m{Zra,y - I3{y - Zi))\< - 53 ^, and \Vy^^^m{Zm,y - I3{y - Zi))\< - 


47 r n- 


dvr n 


2-/3 


2 a 


where c depends only on n and n{zm)- Combining these estimates with (3.23) of Lemma 3.4 
we derive the inequalities 

c c 

\G^{zm,y - I3{y - zi))\ < -and \VyG^{zm,y - I3{y - zi))\ < - 


dvr n— 


dvr n 


2-/3 


2 a 


2 • 


(3.27) 


Then, 


\Rl{Zmt I/) I ^ 


^ \ 2^ / 

27r na s \na a Jo 


\Ua{y -/3{y - zi))\d/3 + J \VyUa{y -/3{y - zi))\d/3j . 

(3.28) 


Then, for I ^ m, (3.21) and (3.28) and observing that Gi is a constant in imply the 
estimate 


G,{zm, y)K^{y)Caiy)Ua{y)dy - G^{zm, zi)K^{zi)CiUa{zi)Vol{ni) 


< 


cCiKaizi) f 

2(2-13) / 

TT 3 Jni 


\Ua{y - ld{y - zi))\dld \y - zi\dy 


+ 


cClKaizi) f 

2-13 I 

271 na 3 Jui 

< [Ka{zi)]Cl i(2-p) < it'maxC 2(2-/3) 

« ^ 3 , 


\VyUa{y-(3{y-zi))\dl3 


p.iej'-l 2(2-13) 

n^a 3 


\y- zi\dy 
(3.29) 


for a suitable constant Ci. 

Regarding the integral G^izm, y)Ca(y)Ga(y)dy we do the following estimates: 


G^(zm, y)Ka{y)Ca{y)Ua{y)dy 


3 ^ M \ 3 2(2-P) 

< ^^C^Krr^a^C {\ a 3 , 


(3.30) 


From (3.20), we can have 


| 24 |<| / G^{zm,y)Ka{y)Ca{y)Ua{y)dy\ 
J 
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[a/5-2] 

E 

1=1 

j¥=m 


\G,,{Zm,Zj)Ka{Zj)CjUa{Zj)Vol{Qj) 



G^{zm, y)Ka{y)Ca{y)Ua{y)dy\ 


which we can estimate by 


[2a 


/3-2, 


w< E 2[(2n + 1)^ - (2n- 1)3] 


n=l 


\G^{Zm,Z,)Ka{z,)C,Ua{z,)V0l{n,)\ 


\G^{zm,y)Ka{y)Ca{y)Ua{y)\dy 


+ \ G^{zm,y)Ka{y)Ca{y)Ua{y)dy\. 


and then 
Finally 


I ^ r 2(2-/3) 2-13 

|A| < ^ -|-a 3 ]. 


1^1 < GK^axa ■ 


3.4.2 Estimate of B 

M 

G^{z^, z,)C,Ua{z,)a^-^ - 


[Zm, Zj)Ka{Zj)GjUa{Zj)Vol{Qj) 


1=1 


1=1 

j¥=m 


[KaiZm)] 


K-2] [Kaiz,)] 


G^{Zm,Zl)GlUa{zi)a‘^ X] X] G^,{Zm,Zl)GlUa{zi)c 


,2-/3 


^ = 1 
Ij^m 

Zi 


j=l 1=1 

j^m zi^Q,j 




- G^{Zm,Zj)Ka{Zj)GjUa{Zj)Vol{Qj) 


1=1 

iT^m 


[KaiZm)] 

G^{Zm,Zi)Ua{Zi) 


1 = 1 
l^m 

Zi 


[Ka{Zj)\ 


+ ^Gja^'^[{ G^{Zm,Zl)Ua{zi)) - G^{Zm,Zj)[Ka{Zj)]Ua{Zj)\, 

1=1 /=! 
j^m zi^Q.j 

since VoliVtj) = a?~^ Ci = Cj, for I = 1,..., Ka{zj). We write, 


[KaM] 

E{ := G^{Zm,Zi)Ua{Zi) 

1=1 

l^m 

ZlG^m 


(3.31) 
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[Ka{Zj )] 

E 2 • ^ ^ GZl^U [■^ci('^j)]^cc('^j)] 

l=l 

Zi^^j 

[Kaizj)] 

= {GZZm,Zi)Ua{Zi) - GZZm,Zj)UaiZj)). (3.32) 

l=l 

Zi^^j 

We need to estimate Cma^~^E{ and ^ Cja?‘~^E{. 

j^m 


Now by writing f{zm,y) ■= G^Zm, y)Ua{y)■ For zi e Qj, j 7^ m, using Taylor series, we 
can write 

f i^Zmi Zj) f (Zffi, Zij (^Zj z{)R i^ZiYi^ Zj, Zij, 

with 


R {Zfn^ZjfZl') / ^yf ^Zfjii Zj l^{,Zj Zl')')di(3. 


(3.33) 


By doing the computations similar to the ones we have performed in (3.22 3.28) and by using 
Lemma 3.3|and Lemma 3^ we obtain 


C,a^-f^Ei\ < C2CKma.a-^ (3.34) 

i=i 

j^m 


One can easily see that, 


1/^ „2-/3z7Ti| ^ Ci{Kmax 1 )C O ^ _ Ci{Kmax 1 )F ^2-/3-t 

1^"*“ *iis ^ « ■ 


(3.35) 


Substitution of (3.19) in (3.18) and using the estimates (3.29) and (3.30) associated to A 


and the estimates (3.34) and (3.35) associated to B gives us 


M 


Ua{Zm) + J2^-(^rn,Zj)CjUZZj)a^-^ = -V^Zm^O) 


(3.36) 


i=i 

j^m 


+0 (c^CK^^a^) + O j , 


We rewrite the algebraic system (2.9) as 

M 


Ua,m + G^Zmi Zj)GjUa,jO? ^ — —V^Zm) 


(3.37) 


i=i 
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where we set Uam '■= CrnQnn recalling that Cm = C'm o? ^. 


Taking the difference between (3.37) and (3.37) prodnces the algebraic system 

M 

{Ua,m-Ua{Zm)) + Y,G^{Zm,Z^)C,{Ua,j-Ua{Zj))a^-^ = O {CKma.{a^ + 
i=i 


Comparing this system with (2.9) and by using Lemma 3.1, we obtain the estimate 


M 

J2{Ua,m-Ua{zm)) = O 


(3.38) 


m=l 


For the special case d = a^, M = 0{a^ with t > 0, we have the following approximation 


of the far-held from the Foldy-Lax asymptotic expansion (2.6) and from the dehnitions Ua,m ■= 
CmQm and Cm := -\m\dDm\ = Cma^~^, for m = 1,..., M: 


M 


47rU^{x,e) = v:^ + Y,Vn{zj,-x)C,Ua,,a^-^ + 0{a^-^-^f^) 


i=i 

M 


Vn{zj, -x)C,Ua, + O 
j=i 


Consider the far-held of type: 


U^^{x,9) = v:- + — Vr,{y,-x)K,{y)Ca{y)U,{y)dy. 


(3.39) 


corresponding to the scattering problem (3.15). Taking the diherence between (3.40) and (3.39) 
we have: 


A7riU^^{x,e)-U^{x,e)) 

' Vn{y, -x)Ka{y)Ca{y)Ua{y)dy - Vn{zj, -x)C,Uay-^ + O 


M 


>Q 

r7./3-2 


j = l 


K-2] [KM^)] 

I Vn{y,-x)Kaiy)Ca{y)Ua{y)dy- E V^izi,-x)QUa,ia^-f^ + O 


j=i 


iP-2] 


j=i 1=1 

Zl^^j 


V Ka{zj)Cj / [Vn{y, -x)Ua{y) - Vn{zj, -x)Ua{zj)] dy 

j=l 




,2-4 


J=1 


[Ka{zj )] 


[Ka(zj)] 


(VniZj, -x)Ua{Zj) - Vn{zi, -x)Ua{zi)) + K(^Z, -x) {Ua{zi) - Ua,l) 


1 = 1 
ZiGflj 


1=1 

ziGQ,j 
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+ O (ai-3^) 

V / Ka{zj)Cj [Vn{y, -x)Ua{y) - Vn{zj, -x)Ua{zj)] dy 

j=l 

[a/ 3 - 2 ] [Kaizj)] M 

+ ^ {Vn{Zj,-x)Ua{Zj) -Vn{zi,-x)Ua{zi)) + J2yn{Zj,-x)Cja^-^ [Ua{Zj) 

j=l 1=1 j=l 

zi&nj 


Ua,] 


+ 0{a^-^^) 


K-2] 

j3.38| l ^ ^ j I [V^(|/, x'jUai^y') Vni^Zj, x'jUai^Zj')^ dy 

J = 1 
r„/3-2i 


Qj 
lKa{zj)] 




j=i 


1=1 

ZiGQj 


+ O (ai-3^) 


(3.40) 


Now, let us estimate the difference ^ ^a{zj)Cj L [14(2/, ~^)Ua{y) — Vn{zj, —x)Ua{zj)] dy. 

J 

Write, fi{y) = Vn{y, —x)Ua{y). Using Taylor series, we can write 


with 


fi{y) - fi{zj) = {y-Zj)- Rj{y), 


Rjiy) = / ^yifi)iy - Piy - Zj))dp 


[Vy [Vn{-x, y- /3{y- Zj))Ua{y - f3{y - zj))]] d/d 
[VyVni-x, y- f3{y- Zj))] Ua{y - f3{y - Zj)) did 


+ / Vni-x,y - /d{y - Zj))[VyUaiy - ld{y - Zj))] d^. 
Jo 


(3.41) 


Recall that 14 satisfies the scattering problem (2.7)- (2.8), hence it is also solution of the 
corresponding Lippmann-Schwinger equation Vn(x) -|- f^^(x, y)(n^(y) — l)Vn(y)dy = — 


This is the same integral equation as ( 3.12 ) at the expense of replacing KaCa by n — 1 and 14 
by Then replacing in Lemma 3^ 14 by and then Y by 14, we have the following 

estimates 

||14|4 oo(q), II V14||l°°(o) 4 C (3.42) 

where C depends only on ||n||£,00(^7) and n. Then 


\RM\ < c[ nu,{y-/d{y-zy))\d/d + [\vyU,{y-ld{y-z,))\dld]. (3.43) 
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Using (3.43) we get the estimate 


^ Ka{zj)Cj / [Vn{y, -x){y)Ua{y) - Vn{zj, -x)Ua{zj)]dy 




<CJ 2 Ka{zj)CjK / 

j = l h 


\y-^j\ / \Ua{y - f3{y - Zj))\df3dy 


+c Y. K4zi)a 

i=i 




M / \y-^j\ \^yUa{y - ^{y - Zj))\d^dy 

Qi Jo J 


) 


a/3-21 


<C Ka{zj)CjCia^ ^0 3^ _l_ Cg) 

j=i 

/ \ 2-/3 

<CK,naxCCi[K + C5)a 3 . 

In the similar way, using ( |3.38[ ), we have, 

a^-2] [Kaizj)] 

Y, iyn{Zj,-x)Ua{Zj)-Vn{zu-x)Ua{zi)) 
j=l 1=1 

Zi ^ 


(3.44) 


^ O ( K-maxC-i^d 3 Oj 




(3.45) 


Using the estimates (3.44) and (3.45) in (3.40), we obtain 


-u^^{x,e)-u^{x,e) 

= O {x^axa^Cci (« + cs)) + 0(C(C + l)K^axM + a^"^"*)) + O 

= 0(a^ + (3.46) 

Since VoI{VL) is of order + i)^) S'lid d is of the order a*, we should have t > 

Hence, we need to impose the following conditions 

2 — B 

t > -, 2 — f3 — t > 0 and 1 — 3/? > 0 . 


3.5 End of the proof of Theorem 2.2 


Combining the estimates (3.46) and (3.17), we deduce that 


— [U^{x,e)-U^{x,e)]-x = 0{a 


- V- 2-/3-t, 1-3/3})^ ^ ^ < t < 2 -/3 


(3.47) 


uniformly in terms oi x, 6 G §^. 
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